Excluding Light Asymmetric Bosonic Dark Matter 
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We argue that current neutron star observations exclude asymmetric bosonic non-interacting 
dark matter in the range from 2 keV to 16 GeV, including the 5-15 GeV range favored by DAMA 
and CoGeNT. If bosonic WIMPs are composite of fermions, the same limits apply provided the 
compositeness scale is higher than ~ 10 12 GeV (for WIMP mass ~ 1 GeV). In case of repulsive 
self-interactions, we exclude large range of WIMP masses and interaction cross sections which com- 
plements the constraints imposed by observations of the Bullet Cluster. 
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1. Introduction. An appealing solution to the dark 
matter (DM) problem is offered by Weakly Interacting 
Massive Particles (WIMPs) emerging in many theories 
beyond the Standard Model (SM). However, WIMPs are 
very difficult to detect, and therefore little is known 
about their properties. Experimentally, the situation is 
rather unclear (see e.g. limits from CDMS with 
DAMA and CoGeNT suggesting the existence of a 
light WIMP with a mass around ~ 10 GeV. 

Apart from direct searches, constraints on WIMPs can 
be set by observations of compact objects such as white 
dwarfs and neutron stars (j-QjJ. These constraints can 
be grouped in two types. The first type targets WIMPs 
that can annihilate inside the star producing heat that 
can change the thermal evolution of the star [5J . WIMPs 
of this type can arise in supersymmetric extensions of 
the SM (see 1131] and references therein), or in Techni- 
color models [1J|- Constraints of the second type tar- 
get asymmetric DM models. In these models the anni- 
hilation of DM in the present-day Universe is impossi- 
ble because only particles, and no anti-particles (hence 
the term "asymmetric") remain [ToT - fTij . An additional 
bonus in these models is that the asymmetry of WIMPs 
might be linked through sphalerons with the baryon 
asymmetry [l6l Il8j|. which can explain the today's ratio 
Udm/^b ~ 5 provided the WIMP has a mass around 
5 GeV. Note that this value is not far from the one sug- 
gested by DAMA and CoGeNT. In view of this coinci- 
dence, the models with WIMP masses in the GeV range 
have become quite popular. 

Since in the asymmetric DM models WIMPs cannot 
annihilate, if a large number of them is accreted during 
the lifetime of a neutron star, they may collapse forming 
a small black hole inside the star that eventually destroy 
the latter. Therefore the existence of old neutron stars 
can impose constraints on the properties of asymmet- 
ric WIMPs. In fact, in the case of fermionic asymmet- 



ric WIMPs with a spin-dependent cross section, these 
constraints are competitive to direct DM search experi- 
ments [TTj j . 

In this letter we focus on asymmetric bosonic dark 
matter and derive constraints on the DM parame- 
ters from the formation of black holes inside neutron 
stars. We show that for fundamental asymmetric non- 
interacting bosonic WIMPs, current observational data 
exclude all the masses from 2 keV to 16 GeV (includ- 
ing the range of masses 5 — 15 GeV favored by DA VI A 
and CoGeNT). If WIMPs are composite particles made 
of fundamental fermions, the above constraint does not 
apply. However, if the compositeness scale is above 
- 10 12 GeV (like in Grand Unified Theories (GUT)), 
candidates like that are again excluded in the same mass 
range. In addition, we constrain the case of fundamental 
self-interacting bosonic WIMPs. We show that if the in- 
teraction is repulsive, a vast area of self-interaction cross 
sections complementary to the one excluded by observa- 
tions of the Bullet Cluster [2(| is excluded. 

2. Bosonic Dark Matter. Gravitational collapse of a 
self-gravitating lump of particles happens differently in 
case of bosons and fermions. In the case of fermions of 
mass to, a large number of particles N ~ (Mpi/to) 3 is re- 
quired to overcome the Fermi pressure, where Mp\ is the 
Planck mass. In the case of non-interacting bosons this 
number is parametrically smaller, N ~ (2/7r)(Mpi/m) 2 , 
since only the uncertainty principle opposes the collapse. 
A repulsive interaction of bosons would provide an extra 
pressure, so the required number of particles is larger in 
this case. Taking a X(f> 4 model as a generic example, the 
minimum mass of a self-gravitating lump which can form 
a black hole is [H[ 
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where we have expressed the result in terms of the self- 
interaction cross section a = A 2 /(647rm 2 ). Here and be- 
low we use the natural units H = c = fcs = 1. 

It is easy to see from eq. (Q]) that for cross sections 
a ^> Mpjm 2 ~ 10~ 104 cm 2 (m/GeV)~ 2 , the second term 
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dominates, and the minimum required mass scales (at 
constant A) in the same way as for the fermionic parti- 
cles with a different (and potentially much smaller) co- 
efficient. The best experimental constraints on the self- 
interaction cross section come from the Bullet Cluster, 
ct/to < 2 x fO- 24 cm 2 /GeV 0. 

Several conditions have to be satisfied for a gravita- 
tional collapse of WIMPs inside a neutron star to occur. 
Firstly, a sufficient number of DM particles must be ac- 
cumulated during the lifetime of the neutron star. The 
accretion of WIMPs onto a typical 1.4M neutron star in 
a globular cluster, taking into account relativistic effects, 
has been calculated in [9(. The total mass of accreted 
WIMPs is 

M "» = «* 1046 (iod^)fe)' GeV ' <2 > 

where the "efficiency" factor / = 1 if the WIMP- 
nucleon cross section satisfies a n > 10 _45 cm 2 , and / = 
cr (1 /(I(r 45 cm 2 ) if cr„ < ICT 45 cm 2 . The condition 

M acc > M cr u (3) 

guarantees that the accumulated DM mass is above the 
critical value ([!]). 

Secondly, the newly-formed black hole must accrete 
matter faster than it evaporates due to Hawking radia- 
tion. In the Bondi regime of accretion, the change of the 
black hole mass M with time is given by the equation 

dM _ 4tt Pc G 2 M 2 1 

~dT ~ cf I53607rG 2 Af 2 ' ^' 

where c s and p c are the speed of sound and the mass 
density of the neutron star core, respectively. The first 
term corresponds to the Bondi accretion while the sec- 
ond represents the energy loss due to Hawking radiation. 
Since the accretion increases while the Hawking radiation 
decreases as a function of M, it is the initial mass of the 
black hole that determines its fate. Requiring that the 
first term dominates when the black hole is formed gives 
the condition 

M > 5.7 x I0 36 GeV. (5) 

Here we have used p c = 5 x 10 38 GeV/cm 3 and c s = 0.17. 
Any black hole with the initial mass satisfying eq. ([5]) will 
eventually destroy the whole star, while the smaller black 
holes will evaporate with no detectable effect. 

The third condition necessary for the WIMP col- 
lapse into a black hole is the onset of the WIMP self- 
gravitation. WIMPs captured by the neutron star ther- 
malize within a time t th = 2x 10 _5 yr (m/GeV) 2 0,11,11] 
and concentrate in the center within the radius 




where T c is the temperature of the star core. When their 
total mass M increases beyond the mass of the ordinary 



matter within the same radius, 

A f iryi \ — 3/2 

M >- Wt 3 h = 2.2xl0 46 GeV(— ) , (7) 

their own gravitational field starts to dominate over that 
of the star and the self-gravitation regime sets in, leading 
to the gravitational collapse provided the condition ([3]) is 
satisfied. It can be seen from eq. © that ([7]) is satisfied 
if the WIMP mass is larger than -I GeV (-143 GeV) 
for p dm = I0 3 GcV/cm 3 (p dm = 0.3GeV/cm 3 ), but not 
for lighter WIMPs. 

However, if WIMPs are bosons they can form a Bose- 
Einstein condensate (BEC). Since this state is more com- 
pact, the self-gravitation in this case starts for a smaller 
number of particles, i.e., before the condition (J7]) is sat- 
isfied. The particle density required to form BEC is 

,„ o / m \3/2 / T c \ 3/2 

n ~ 4.7 x 10 cm ( — irr 

\GcVJ \W 5 KJ 

Assuming an old neutron star with a temperature T c — 
10 5 K, the number of WIMPs needed in order for BEC 
to form is iV BE c - 2 X 10 36 . All the WIMPs accreted in 
excess of this value will go into the condensed state. For 
most of the cases of our interest, the number of accreted 
WIMPs will be larger than Vbec, so eq. ([7]) has to be 
reconsidered. 

The size of the condensed state is determined by the 
radius of the wave function of the WIMP ground state in 
the gravitational potential of the star, 

/8tt 2 V 1/4 4 /GeV \ 1/2 
r c = ( —Gp c m 2 J ~ 1.6 x 10~ 4 ( J cm. (8) 

Substituting this size in place of rth in eq. ((7J) we get 

M> 8 x 10 27 Gev(^)~ 3/2 . (9) 

In view of eq. ([2]), the amount of DM sufficient for WIMP 
self-gravitation in the condensed state can always be 
accumulated provided that the WIMP is heavier than 
— 0.1 eV, which covers all cases of interest. Thus, due to 
the formation of BEC the requirement of self-gravitation 
does not provide an extra condition. 

Finally, the accumulation of WIMPs may become in- 
efficient if they may escape from the neutron star once 
captured, which is a danger at small WIMP masses. It 
can be seen from eq. © that for WIMP masses in the 
keV range the radius of the WIMP lump becomes com- 
parable to the size of the star, so that WIMPs in the tail 
of the velocity distribution may escape. The rate F of 
WIMP evaporation can be estimated as follows [23 |. 

F = n s f^V'Vl + ^W-^V 
\2irm ) \ RT ) v \ RT J 

(10) 

where T, M and R are the temperature, mass and radius 
of the star, respectively, and n s is the WIMP density at 
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FIG. 1: The exclusion regions as a function of the DM mass 
and the WIMP-to-nucleon cross section for an isolated neu- 
tron star at local DM density p = 0.3GeV/cm 3 (such as 
J0437-4715 and J0108-1431) and for a neutron star in the 
core of a globular cluster at p = 10 3 GeV/cm 3 . 

the surface. Calculating n s from the Boltzmann distribu- 
tion we found that the evaporation can be safely ignored 
for masses larger than ~ 2 keV. 

In summary, accumulation and subsequent gravita- 
tional collapse of WIMPs captured inside a neutron star 
occur for WIMPs heavier than ~ 2 keV provided the 
conditions (|3J and © are satisfied. In the case of no 
self-interactions the collapse to a black hole inside the 
neutron star happens for WIMP masses 2 keV < m < 
16 GeV. The upper bound of this exclusion range is in- 
dependent of the local DM density p, while the lower 
bound raises slightly at small p. As a function of the 
DM mass, the bound on the WIMP-to-nucleon cross 
section is a n > 8 x 10 _50 cm 2 (GeV/m) for nearby iso- 
lated stars at local DM density p — 0.3 GeV/cm 3 and 
rj„ > 2 x 10 _54 cm 2 (GeV/m) for stars in globular clus- 
ters at local DM density p = 10 3 GeV/cm 3 . The resulting 
exclusion regions are shown in Fig. [TJ 

Several old neutron stars have been observed, both in 
the vicinity of the Earth where the DM density is p ~ 
0.3 GeV/cm 3 , and in the cores of globular clusters where 
the DM density may be as high as 10 3 — 10 4 GeV/cm 3 
Q. Examples of nearby neutron stars are J0437-4715 [24| 
and J0108-1431 [l| (140 pc and 130 pc from the Earth, 
respectively). The examples of neutron stars in globu- 
lar clusters are e.g. the pulsar B1620-26 located at the 
outskirts of the core of M4, and X7 from 47 Tuc [Hj]. 

In the case of a repulsive interaction, the exclusion re- 
gion that follows from eqs. (j3|) and (0 is shown in Fig. [2] 
Depending on the self-interaction cross section, the con- 
straints extend to much higher masses and are comple- 
mentary to those derived from the observation of the Bul- 
let Cluster. 

3. Composite Dark Matter. The discussion above 
refers specifically to fundamental bosonic DM. If instead 
the latter is composite of fermions, the situation might 
change. There are two possibilities. Assume the total 
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FIG. 2: Constraints on bosonic DM mass m and self- 
interaction cross section a from observations of neutron stars 
in globular clusters. Excluded region (pink) is shown for two 
background DM densities as indicated on the plot. The cyan 
region shows constraints from the Bullet Cluster. 



WIMP mass exceeds M crit , eq. (JTJ). As the DM lump 
shrinks towards its Schwarzschild radius, it might reach 
the density at which the mean distance between WIMPs 
is comparable to the scale of compositeness. At this point 
the Fermi pressure comes into play and might stop fur- 
ther collapse unless the DM lump has already reached its 
Schwarzschild radius. 

To estimate the minimum compositeness scale A cr i t , 
we express the mean distance d between WIMPs in 
the nearly-collapsing (i.e., having size comparable to its 
Schwarzschild radius) DM lump in terms of its mass 
M. Ignoring the numerical coefficients, we have d — 
GM^/^m 1 ^ . Taking the mass to be equal to the crit- 
ical one, eq. (pTJ) , we get 

A - = ™ V3 < 3 ( 1 + 5 S)" /3 - (id 

In the non-interacting case this gives A cr n — 2 x 
10 12 GeV^/GeV) 1 / 3 , which is well below the GUT mass 
scale of order ~ 10 16 GeV for all masses of interest (cf. 
Fig. [2]). Thus, our constraints are also valid for com- 
posite WIMPs with the compositeness scale higher than 
-10 12 GeV. 

4- Discussion and conclusions. Two remarks are in 
order. In the above analysis we have assumed that the 
black hole that is formed inside a neutron star and is 
not destroyed by the Hawking radiation eventually con- 
sumes the whole star. However plausible, this assump- 
tion requires some justification. In eq. (U) we have taken 
Bondi accretion, which is very efficient and would indeed 
consume the whole star much faster than in 1 Gyr (see 
Ref. [ll[ for the calculation). However, the accretion may 
be slowed down by the angular momentum of the star. 
The Bondi regime cannot be maintained if the angular 
momentum of the falling matter exceeds the one it would 
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have on the innermost stable orbit. One can show that, 
in the absence of momentum transfer, for a typical pe- 
riod of an old neutron star of the order of a second, only 
a small inner part of the star can be consumed in the 
Bondi regime. The momentum transfer has been studied 
in Ref. 26] in the context of ordinary stars. Rescaling the 
parameters to the case of a neutron star, we found that 
if the momentum transfer due to the viscosity is taken 
into account, the Bondi accretion is maintained until the 
black hole reaches a mass of ~ 1O~ 7 M . From this state, 
with the Bondi rate the consumption of the whole star 
would take about 1 min. Even if the actual rate is many 
orders of magnitude slower, the star will definitely be 
destroyed within 1 Gyr. 

The second remark concerns the validity of the Bondi 
regime at the initial stages of the black hole growth. The 
lightest black hole relevant for our analysis has size and 
Hawking temperature in the GeV range. While the fluid 
approximation should still be adequate at these scales 
because of the extremely high density of the nuclear mat- 
ter, the actual parameters may differ from those used in 
cq. ([I]). This might change slightly the upper value of the 
exclusion mass range. The precise calculation is difficult 



and goes beyond the scope of this letter. 

To conclude, we have demonstrated that the exist- 
ing observations of old neutron stars, both in globu- 
lar clusters and in the vicinity of the Earth, exclude 
light non-interacting fundamental bosonic dark matter 
candidates in the mass range from 2 kcV to about 
16 GeV with WIMP-to-nucleon cross section a n > 8 X 
10 _50 cm 2 (GeV/TO). The constraints equally apply to 
composite bosonic dark matter if the compositeness scale 
is higher than ~ 10 12 GeV. A wide range of masses 
and cross sections is also excluded for very weakly self- 
interacting bosonic candidates. 

The work of P.T. is supported by the IISN project 
No. 4.4509.10 and by the ARC project "Beyond Einstein: 
fundamental aspects of gravitational interactions" . 

Note added: When this paper was being finalized 
Ref. [13] appeared which addressed the same question. 
However, in Ref. 27} the Hawking evaporation of black 
holes was disregarded which invalidates the constraints 
for WIMP masses m > 16 GeV. At smaller masses our 
results are somewhat different because the effect of the 
Pauli blocking was overestimated in Ref. [27j |. 
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